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Abstract. In this paper, we are interested in the free Jacobi process starting 
at the unit of the compressed probability space where it takes values and 
associated with the parameter values A = 1,6 = 1/2. Firstly, we derive a time- 
dependent recurrence equation for the moments of the process (valid for any 
starting point and all parameter values) . Secondly, we transform this equation 
to a nonlinear partial differential one for the moment generating function that 
we solve when A = 1,6 = 1/2. The obtained solution together with tricky 
computations lead to an explicit expression of the moments which shows that 
the free Jacobi process is distributed at any time t as (l/4)(2-|- Yijt +^2*) where 
Y is a free unitary Brownian motion. This expression is recovered relying 
on enumeration techniques after proving that if a is a symmetric Bernoulli 
random variable which is free from {Y, Y*}, then the distributions of Y2t and 
that of aYtaY* coincide. We close the exposition by investigating the spectral 
distribution associated with general sets of parameter values: A = 1, S G (0, 1) 
and A G (0,2), 6 = 1/2. 



1. Reminder and motivation 

The one dimensional Jacobi process is a two parameters-dependent stationary 
random motion valued in the interval [0, 1] and its equilibrium measure is given 
by the Beta distribution ([8]). For special values of the parameters, this process 
can be realized as norms of projections of a spherical Brownian motion in a finite- 
dimensional Euclidean space onto spheres of lower dimensions {[1^). When extended 
to matrix spaces, this realization gives rise to matrix Jacobi processes as radial 
parts of corners of Brownian motions in the orthogonal and the unitary groups, 
and the corresponding equilibrium measure is the matrix- variate Beta distribution 
([Z], [in])- Motivated by free probability theory, the free Jacobi process was realized 
in [9] as the rescaled limit (in the sense of noncommutative moments) of the complex 
matrix Jacobi process. Its equilibrium distribution was determined in [6^, and 
[5] . By the virtue of asymptotic freeness of independent unitarily-invariant random 
matrices and constant ones whose distributions converge, the matrix model of the 
free Jacobi process suggests the following abstract definition (see [9] for details). 
Consider a non commutative W*-probability space (^,t), that is a von Neumann 
algebra with unit 1 and endowed with a faithful normalized trace r. Let 9 G 
(0, 1), A > such that < X9 < 1 and take two projections P, Q in jz/ such that 

• T{P) = X9,T{Q)=e, 

• PQ = QP = P {P < Q) a X < 1 and PQ = QP = Q {P > Q) otherwise. 
Let also K be a free unitary Brownian motion in ^ ([3]) and assume 



{p,g}, {Y,Y*} 
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are free families in si . Then 

Jt := PYtQYtP 

defines the free Jacobi process of parameters (A, 0) and one easily sees from 

that P — J is a free Jacobi process too of parameters (A6'/(l — 6*), 1 — 0). However, 
the stochastic analysis of J performed in [9 requires that Jo and P — Jo are injective 
operators in the compressed probability space 

For that reason, the free Jacobi process studied in [9] and denoted J as well is driven 
by a free unitary Brownian motion starting at Z, where Z is a unitary operator in 
£/ such that 

{P,Q},{Z, Z*},{Y,Y*} 
are free families in s/. Actually, the operator Z has to ensure that 
T := inf{i > 0, Jt orP — Jt is non injective} 

is positive and the issue of our stochastic analysis is the free stochastic differential 
equation 

(1) 

djt = Vxe^/7tdWt\/p - Jt + Vx0^/7tdWt y/p~~Jt + {ep - Jt)dt, o < t < t, 

where 14^ is a Pjz/P-complex free Brownian motion. This equation was the key in- 
gredient for investigating spectral properties of Jt . Namely, let TO„ (i) := t{JI^)/t{P) 
be the moments of Jt in Ps/P then for any n > 1 and up to time T ([9]) 

n-2 

(2) dtmn{t) = -nrrinit) +6'nm„_i(<) + X9n^^mn-k-i{t){mk(t) - mk+i{t)) 

where the sum in the RHS is taken to be empty when n ~ I and mg (t) ~ I, nin (0) = 
t{{PZQZ*)^)/t{P). Equation may be thought of as a weak version of ([1]) and 
since Jt exists at any time t and for any starting point Jq then it is very likely that 
^ remains valid for any initial data m„(0),n > 1 (that is independently from Z). 
In this paper, we shall confirm the validity of ([2|) relying on the time-dependent 
recurrence equation derived in ^2 , Theorem 3.4 for 

T{aiYta2Yt . . . a2n~iYta2nYt), 

where {oi, . . . , a2n} is a 27i-tuple of random variables forming a ^-free family with 
Y . More precisely, we first specialize this theorem to a/c = P when k is odd and 
a-k = Q when k is even, then infer that the computations keep unchanged when we 
rather consider Uk = Z*PZ when k is odd and Z is any unitary operator in which 
is *-free from Y and from {P, Q}. Indeed, r is a trace, Z*PZ and P have the same 
moments and both families {Z*PZ, Q} and {Y, Y*} are free. Doing so gives us the 
benefit of studying the free Jacobi process PYQY*P which was discarded in [5] and 
which is related to Voiculescu's liberation process associated with free projections 
f[16|). However, observe that for any A, the sequence Vn{t) := Xmn{t),n > 1 
satisfies 

n-2 

dtVn{t) = -nVn{t) +enVn-l{t) + 9n^Vn-k-lit)iVkit) - Vk+l{t)) 

k=a 
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with initial data Vn{0) — Am„(0), vo{t) = Amo(i). Thus, we shah primarily consider 
the value A = 1 and focus on the corresponding free Jacobi process. However, due 
to the high nonlinearity of ©, we further restrict our attention to the value 6 = 1/2 
and prove that the spectral distribution of Jt in the compressed space, say fit-, fits 
the distribution of the random variable 

l[Y2t+Y,\ + 21] 

in (^,t). In particular is absolutely continuous with respect to Lcbesgue mea- 
sure on the real line and its support fill in the interval (0,1) at time t = 2. Actually, 
our description follows from a closed formula for m„(t),n > 1, namely: 

where is the n-th Laguerre polynomial of index 1 ([15 ). Formula ([3]) is in turn 
obtained after solving a nonlinear partial differential equation (hereafter p.d.e) for 
the moment generating function of /i^; 

Mt(z) = ^TO„(t)z" |z|<l. 

n>0 

In fact, by the virtue of this p.d.e admits a unique solution in the class of 
analytic functions around zero. Nevertheless, the binomial coefficients present in 
(jS]) suggests that this formula might be derived based on enumeration techniques. 
To this aim, we start by noting that when A = 1, the faithfulness of r and P < Q 
entail P = Q: indeed t(P) = t{Q) so that 

r{{P - Qf) = 2t{P) - 2t{PQ) = 0. 

Next, if further ~ \/2 then P = (l + a)/2 where a € is a self-adjoint symmetric 
Bernoulli random variable and we shall derive the following expansion: 

2r((l + a)Y,{l + a)r/)") = (^^j) + ^^"^ ^r{{aY,aY:f). 

By the way, we are grateful to A. Nica and to C. Sattlecker for their help in 
completing the proof of this derivation. Finally, we recover ([3]) after proving that 
Y2t and vt '■= oYtoY^ have the same distribution. The enumeration techniques not 
only provides another elegant proof of the description of /i*, but opens the way 
to investigate the spectral distribution of Jt for arbitrary 9 S (0, 1) and A = 1. 
More precisely, the above expansion of r((l + a)Yt{l + a)!^^*)") comes in with the 
additional term 2^"~^(20— 1). However one has to describe the spectral distribution 
of aYtoY^ which does not seem to be easy. As to the parameter values A € (0, 2), 6* = 
1/2, we already noticed above that they correspond to a deformation of the initial 
data TO„(0), mo(t). However we should consider the parameter values A G (0, 1] and 
A € [1,2) for fixed 9 = 1/2 separately, since PQ = P ^ to„(0) = 1 and PQ = 
Q mn{0) = 1/A respectively. Nonetheless, note that the moments (m„(f))n 
corresponding to A € [1,2) can be recovered from the ones corresponding to A € 
(0, 1]. Indeed, if P is a projection of continuous rank t{P) = A/2 then 1 — P is a 
projection of continuous rank (2 — A)/2. Moreover, since YtQY^ is idempotent and 
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since r is tracial then we easily derive 

r[((l - P)Y,QYm = r(Q) + ^(-1)'= ( ^ ) r[(Pr,Qr,*)^-] 

k=l ^ ^ 

Accordingly, we only need to consider the parameter values Ae(0,l],6' = l/2 and 
we shall prove that the moment generating function of the spectral distribution 
of Jt splits into two parts. The first 'principal' part is the sum of the moment 
generating function of the stationary distribution (i.e. the spectral distribution 
in the compressed space of PUQU*P where J7 is a Haar unitary variable in jz/ 
which is T^r-free from {P, Q}) and of the generating function of a linear combination 
of the moments of i^2At/(2-A)- The second part is an analytic function ut around 
the origin whose first and second derivatives vanish there. Unfortunately, we do 
not succeed to provide a description of the spectral distribution of the free Jacobi 
process through known random variables as we already did forA = l,6'==l/2. 

2. A RECURRENCE TIME-DEPENDENT EQUATION 

In this section, we prove that 
Proposition 1. The moments of the free Jacobi process 

mnit) := ^TiPZYtQY^Z^P)"), n > 1 

satisfy ^ for any unitary Z that is -k-free from Y and from {P, Q}. 

Proof. As explained in the introductory part, computations are given in the 
special case Z = 1. Moreover, we shall equivalently prove that the sequence 

rn{t) t(J,") = T{P)mn{t),n> 1 

satisfies 

Tl-2 

(4) dtrn{t) = ~nrn{t) + n9r.n^i{t) + n ^ r„_fe_i(t)(rfe(^) - rk+i{t)) 

fc=0 

with ro{t) = t{P) = \9. To proceed, we recall Theorem 3.4 in [2]: 
Theorem 1. Let n>l and define 

/2„(ai, . . .,a2„,t) := e"*r(aiyfa2F* . . . a2„-iYta2nY*) 
where {ai, . . . , 02,1} € is -k-free with Y . Set faiA, t) := t{A) for any A € then 

dtf2n{ai, . . . , a2n, t) = — ^ /2n-((-fc)("l' • ■ • 7 ^fc, Ol+i, . . . , a2n, t)fi^k{ak+i, ■ ■ ■ ,ai,t) + 

l<k<l<2n 
l-k=a[2] 



n— 

(;_fe)_i(ai,. . .,afc_i,afcai+i,a;+2, . . . ,a2n,t)fi-k-i{aiak+i,ak+2, ■ ■ . ,a;_i,i) 



l<fe<i<2ri 
l-k-l=0[2] 

Now, we specialize Theorem [T] to = P if fc is odd and ak = Q otherwise. 
Then the result is straightforward when n — I and is even stated in [2j, p. 923. 
So let n > 2 and note that both indices k, I in the first (respectively second) sum 
in Theorem [T] have the same (respectively different) parity, therefore k and I + 1 
in the second sum have the same parity and so do I and fc + 1. Accordingly, the 
first sum does not contain terms fo{-,t) while the second does: they correspond 
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to indices I — 2n, k — 1 and to I = fc + 1,1 < k < 2n — I. Since P and Q are 
idempotent and since r is a trace, then the contribution of indices k = 1, / = 2n is 
T{P)f2n-2{P, Q,...,P,Q) while that of / = fc + 1, 1 < fc < 2n - 1 is 

[nT{Q) + (n - l)T{P)]f2n-2{P. Q,...,P,Q,t) 

respectively. Thus, both contributions sum up to 

(5) n(r(Q) + t(P))/2„-2(P, Q,...,P,Q, t). 

Next we write I = k + 2s + I for integer positive values of s and distinguish n = 2 
and n > 3. li n — 2 then there is no additional term in the second sum, while if 
n > 3 we separate k — 1 and 2 < fc < 2n — 3. By the same properties of P,Q,t 
mentioned above, the contribution of indices fc = 1, / = 2s + 2 is 

n-2 

(6) h(n-s-i) {P.Q,---, P, Q, t)f2s{P, Q,...,P,Q,t). 

s=l 

For the remaining values of 2 < fc < 2n — 3, we distinguish even and odd ones: the 
contribution of indices fc 2j. 1 < j < n — 2, ^ = 2j + 2s + 1 is 

n—2 n—j—1 

E f2in-s~l){P,Q,---,P,Q,t)f2siP,Q,...,P,Q,t) 

j=l s=l 

while that of fc = 2j + 1, 1 < j < n - 2, Z = 2s + 2j + 2 is also 

n—2 n—j—1 

Y E f2{n-s-l){P,Q,---,P,Q,t)f2siP,Q,...,P,Q,t). 

j=l s=l 

By rearranging the terms in both obtained double sums, we see that the contribution 
of indices 2 < fc < 2n - 3, ^ = 2fc + 2s + 1 is 

n-2 

2 ^{n - s - l)/2(„-.-i)(P, g, . . . , P, Q, t)f2s{P, Q,...,P,Q,t) 
which simplify after the index change si->n — s— Ito 

n-2 

(7) - 1) E /2(n-s-l) {P,Q,--., P, Q, t)f2s{P, Q,...,P,Q,t). 

Similarly we consider the first sum in Theorem [T] and write ^ = fc + 2s,l<fc< 
2n - 2, 1 < s < n - [(fc + l)/2], n > 2. Then it contributes to 

n-l 

(8) n J2 f2in-s){P, Q,---,P,Q, t)f2s{P, Q,...,P,Q,t). 

s=l 

Now we recall that f2n{P,Q, ■ ■ ■ , P,Q,t) — e'^'^rn{t) and take into account the 
exponential factor in front of the second sum of Theorem [T] If n > 3 then ^ and 
dH]) sum up to 

n-2 n-2 

e"*[-nr„_i(t)ri(i) - ^ r„_,_i(i)r,(i) + n ^ r,(i)[r„_,_i(t) - r„_,(t)]. 

s=l s=l 
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With regard to (O, ([6]), the whole contribution of the RHS of Theorem [T] is 

n-l 

e"*{nr(Q) + n[r(P) - ro(t)]r„_i(0 + n ^ r,(i)[r„_,_i(t) - r„_,(t)]}. 

But dtf2n{t) = e'^-^^ldtTnt + nr„(i)] together with ro(t) = t(P) show that (U) holds 
for any t > and any n > 3. If n = 2 then the whole contribution is given by ([5]) 
and jl]) leading to dU as well. ■ 

3. The case = 1/2, A = 1 

This section is devoted to the description of the spectral distribution /it of Jt 
when A = 1,0= 1/2 and Jq = P. A major step towards it is the following result: 

Proposition 2. Let L^, be the k-th Laguerre polynomial and let 

fe=i 

be the unique solution of f|3], p!4] / 

5tPt + ^S^pt = 0, po{z) = 

m t/ie class of analytic function around zero. Then the moment generating function 
Mt of /it is given by 

"'W = 7r^{' + ^«'((7T=TIp)}'W^'- 

Proof. Before coming through computations, we point out that 

z 

a : z H> 



(1 + VT^)2 

maps the open unit disc into itself. Indeed, the following expansion holds ([ 
p.70) 



1 _ l^(l)„(3/2)„_,„ 



{1 + VT^zr (3)„n! 



Ul < 1 



and is still convergent for z = 1 by Gauss Theorem ([12, p. 49). Hence the expres- 
sion of Mt{z) given in the proposition makes sense for all |z| < 1. Now, let \z\ > 1 
then similar computations leading to Proposition 7.1. in [5] shows that 

G.iz):^'-±^, 
z ^-^ z" 

satisfies the p.d.e. for all (A, &) 

dtGt = {[(1 - 2X9)z - 6(1 - X)]Gt + Xez{z - 1)G?} 

with initial value Gq{z) — l/{z — 1). This p.d.e. simplifies when one substitutes 
A = 1,61 = 1/2 to 

dtGt = {z{z - l)Gl] , Go{z) = 
or equivalently for |z| < 1 

dtMt = -^-d,{{l-z)M^),Mo{z) ^ 



2^^^ '^J' l-z 
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Set Stiz) := VI - zMt{z) - 1, then 

1 



dtSt + z^/Y~zd,St + ^-z^fY~zd,Sl = 



with the initial data 



So{z) = - 1. 

\/l - z 



Next note that a is invertible in a neighborhood of zero with inverse function given 

by 

-1/ N 4z 

(1 + z)-^ 

Keeping in mind that |a(z)| < |z| < 1 in the open unit disc, then a extends to a 
biholomorphic map from the open unit disc onto its image. Moreover, 



z\/l — za'{z) — a{z). 
Hence Ft : z ^ St{a^^{z)) and wt z ^ Ft{e^z) satisfy 

dtF,^zd.F,^\zd.F?^, 



and 

1 
2 

FinaUy one easily checks from 



dtWt + —zdzwf — 0. 



1 + a{z) = ^ 1 — a{z) 



that 

2a{z) 

w^{a{z)) = So{z) = — = 2po{a{z)) <^ wo{z) = 2po{z). 

1 — a(z) 

Since 2p2t and wt satisfy the same partial differential equation as pt , the proposition 
is proved. I 

The moments mn{t) are given by 

Corollary 1. For any n> \ and any t>Q 

Proof: First of all the generalized binomial Theorem yields 

1 _ ^ (l/2)„ 



^ n— 



E^-"' N<i 



where (l/2)n = V{n + l/2)/r(l/2) is the Pochhammer symbol. But Legendre 
duplication formula ([H]) shows that 

(l/2)„ ^ J_ (2n 
n\ 22" 
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Now let 2F1 be the Gauss hypergeometric function then for any k> 1 (^15] p. 70) 

22fe 1 / ^ 

2Fi{ k + -,k + l,2k + l,z 



(fc + l/2)„(fc + l)„ z" 
Using Legendre dupUcation formula again, we derive 

22n+2fe 

(2n + 2k)\ = ^r(?i + fc + l/2)r(n + A: + 1) 

which yields 

/2n + 2/c\ z" _ 22'=r(fc + l/2)r(/c + 1) ^ (A; + l/2)„(fc + 1)„ z" 
n V^r(2fc + 1) (2fc + i)„ ^ 



As a result 



CX3 ^ 00 



2n + 2A:\ 2"+^= 



fe=l n=0 

1 .1 / 2n \ z'^ 



22n+2fc 

fc = l Tt = fe ^ 

The Corollary is proved. ■ 
We are now ready to give the description of [it- 

Corollary 2. When — 1/2, A = 1, the free Jacobi process starting at P is dis- 
tributed in the compressed probability space (Px/ P, 2r) as the random variable 

in is/. Consequently, fit is absolutely continuous with respect to Lebesgue measure 
on the real line and its density is given by 

C 2^arccos(VK)^ 

2M^==^l[o,i](x)dx. 

where kt is the density of the spectral distribution of the free unitary Brownian 
motion. In particular, the support of fit is the whole interval [0, 1] for any time 
t > 2. 

Proof. Let hn{t) — T{Yt"),n > 1 be the moments of Yt in (^,r), then ([3 ) 

g-nt/2 

hnit) = L,\_i(?ii), 



ON THE SPECTRAL DISTRIBUTION OF THE FREE JACOBI PROCESS 



9 



But since Y and Y* — Y ^ have the same distribution then /i„ — h-n, n e Z. As 
a result 

^2n> 



i{C„>'i:(r.)M->-i:(„T.)M-)} 

\ A. — X Ki — X y 



2n 



02n J2 io„ _ i.)^(^2^t ") 



fe=0 ^ 

^^riY,r\i+Y,,rn- 

Finally, the support of /it is entirely determined by the one of fet (see [5]). 
Remark. When \ — 1,0 — 1/2 ([2]) tafces </ie /orm 

n-l 

dtmn{t) + -mn{t) = 2 X! '"^n-i~k{t)[mk{t) - mk+i{t)]. 

Thus since the moments of Yt converge as too to those of a Haar unitary random 
variable, then the moments 

1 /2n^ 
2"2 

satisfy 

n-l 

(9) m„ = ^ m„_fc_i[TOfe - nik+i]. 

f2k\ 

But if Ck ■= (l/(fc + 1)) I is the k-th Catalan number i^|13| ) then Legendre 

\k J 

duplication formula entails 

1 r(fc + i/2) 1 

a result 

n— 1 2 ^i — 1 ^ -j^ n— 1 

mn-k~i[mk - J^fe+i] = ^ ^ k)Cn-k-iCk = -^5^ X! ^"-^-i^fc- 

fc=0 /£=0 fc=0 

Consequently, Q is nothing else but the recurrence relation for Catalan numbers: 

n-l 



Cn — ^ ] Cn-l-kCk- 



k=t) 
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4. Enumerative derivation of the moments 

Recall that if A = 1 then the faithfulness of r forces P ~ Q. If further 9 = 1/2 
then P := (1 + a)/2 where a = a* G is/ is distributed according to 

In the sequel, we shall derive ^ relying on enumeration techniques. More precisely, 

Proposition 3. Let a he a self-adjoint random variable in si distributed according 
to 

and assume a and Y are -k-free. Then for any n > I 

r[{{l + a)Y{l + a)Y;r] = i + E -((«^*«^**)')- 

Proof let 6 := YtaY^* then 

(1 + a)Yt{l + a)Y; = (1 + a)(l + b) 

therefore [(1 + a) (1 + 6)]" consists of words formed by letters picked in the alphabet 
{a, 6} subject to the cancellations {a^'"' = 6^*^ ~ 1, 0^*^+^ = a, fe^^+i = ^j.. Moreover, 
we claim that those formed by an odd number of letters have zero expectation. 
Indeed, any such word has a reduced (after taking into account of cancellations and 
that there is no relation between a and b) expression either aba ■ ■ ■ aba or bab ■ ■ ■ bab. 
The claim then follows from the trace property of r and from T(a) = t(6) =0. As a 
result, we need to enumerate for each n > 1 words {{ab)'', 1 < k < n, (ba)'' , I < k < 
n—1} (see below for the constant term fc = 0). To this end, let c(n, fc), d{n, fc), e{n, fc) 
be the number of words (ab)'' , {ab)^a, [ba)^ respectively. Then the expansion 

[(1 + a)(l + 6)]" = [(1 + a)(l + &)]"-i(l + a + b + ab) 

and the observation 

{abf = [{abf]l = [{abfa]a = [{abf-^a]b = [{abf-\ab) 

give the recurrence relation 

c(n, fc) = c{n - 1, fc) + d{n - 1, fc) + d{n - 1, fc - 1) + c(n - 1, fc - 1), n > 2, fc > 1 
and 

c(i,fc) = 4o + 4i 

where 5ij is the Kronecker symbol. In a similar fashion, 

{abf a = [{abfa\l = [{abf]a = [{abf+^]b = [{abf+^a\{ab) 
gives the recurrence relation 

d(n, fc) = d{n - 1, fc) + c(n - 1, fc) + c(n - 1, fc + 1) + d(n - 1, fc + 1), n > 2, fc > 0, 
with 

d(l,fc) = 4-0, 

while 

[(1 + a)(l + 6)]" = (1 + a + 6 + a6)[(l + a)(l + 5)]"-i 

together with 

{baf = l[baf = a[{abfa\ = b[{abf-^a\ = ab[{baf+^\ 
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give the recurrence relation 

e(n, k) ~ e{n — 1, fc) + d{n ~ l,k) + d{n — l,k — I) + e{n — 1, fc + 1), n > 3, fc > 1, 
with 

e(2, k) = 34o + hi, e(l, fc) = 4o- 

Now, we compute c(n, 0) which correspond to the constant term (not depending on 
time t) in the expansion 

r[((l+a)(l + 6)r]. 

To proceed, we argue that up to a factor, c{n, 0) is the 7i-th moment m„((X)) of the 
stationary distribution of the free Jacobi process 

Indeed, for fixed n > 1 

hm T\{aYtaY*Y\ = T{(aUaU*Y\ 

>-oo 

where J7 is a Haar unitary distributed random variable. But a and UaU* are free 
in and since T(a) = T{UaU*) = then the very definition of freeness implies 
that T[{aUaU*)"'] — for any n > 1. We can also compute c{n,0) by considering 
the von Neumann algebra generated by {a, b} endowed with the state that assigns 
the value 1 to the unit and vanishes otherwise. More precisely, this algebra is the 
free product of the von Neumann algebras generated by {a} and by {b} since the 
latters may be realized as von Neumann algebras of the cyclic group of order two 
so that Theorem 1.6.3. in Tlj applies. Therefore a and b are free there and it 
is obvious that they are symmetric Bernoulli random variable with respect to the 
given state. It follows that c(rt, 0) is the constant term in the n-th moment of the 
two-fold free multiplicative convolution 

^{60 + 62) ^^{60 + 62) 

which may be computed using the S'-transform ([l3]). More precisely, the S- 
transform of this free multiplicative convolution reads 

(^ + 1)^ 
(l + 2z)2 

near the origin. Consequently, the corresponding inverse moment generating is 
given by 

+ 1) _ 1 
(l + 2z)2 ~ 4 

whence we deduce the inverse moment generating of the free multiplicative convo- 
lution 

n=l ^ ' 

Based on the initial values c(n, 0), ci(l, 0), we proceed by mutual induction (on 
k for fixed n then on n) and use the elementary identity 

(n\ ( n \ /n + l\ 



(l + 2z)2 



1 



- 1 
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in order to prove that 

^2n- 1 



c{n, k) = 



n 



- k 



( 2n — 1 \ 
d{n,k) = L_y^„2) 

Note by passing that d{n, k — 1) — c{n, k) which agrees with the recurrence relations 
for d(n, k) and e(n, k). Finally 



c(n, k) + e(n, k) = 



'2n-l\ / 2?i - 1 \ / 2n 
n — k J \n — k — 1 J \n — k / 

In order to recover ([3|), it suffices to observe that a 

Lemma 1. The unitary random variable oYtoY^ is distributed as Y2t- 

Proof, since the distribution of a unitary random variable is entirely determined 
by its moments, we shall prove that the sequence defined by 

Sn{t) ■.^e'''TiiaYtaY;r),n>l, 

is the (unique) solution of 

n-1 

dtSn{t) = -71^ s„_fe(t)sfe(<), n>2, si{t) = 1, 

k=l 

then infer from [14 (Theorem 4 p. 669) that 

Sn{t) - -Ll_,{2nt). 
n 

To proceed, we specialize Theorem [1] to aj. = a for all k and notice that the second 
sum vanishes since = 1 and since r(a) = 0. The contribution of the first sum is 
easily seen to be 

n-1 

-n^ s„_fe(t)sfe(t) 

k=l 

and the value of si(t) = 1 is readily derived from [2] p. 923 using T(a) = 0, T[a?) = 
1. The lemma is proved. I 



5. Further developments: general parameter values 

In this section, we investigate the spectral distribution of Jt for more general 
sets of parameter values: A = 1, e (0, 1), A G (0, 2), 9 — 1/2. 

5.1. Parameter values A = 1,6' G (0, 1). As one easily realizes, the enumeration 
techniques used to expand 

r[((l + «)(l + 6))"] 

remain valid when A = 1,6* e (0,1). Indeed, P = (1 + a)/2 where the spectral 
distribution of a is 

951 + {1-9)5-1, 
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and we need to take into account the contribution of words with odd number of 
letters. By the trace property of r and the relations — ~ 1, this contribution 
is T(a) = {26 — 1) up to a positive integer say c(n) therefore 

r[((l + a)Y,il + a)Y;r] = ^ (^^^ + ^((«^*«^**)') + " 1)^^- 

The coefficient c{n) is easily computed when 9 £ (0, 1) \ {1/2} by letting t = and 
using (1 + a)^ = 4P, 1: 

4"6l = 2^"-! + {26 - l)c(n) ^ c{n) = 2^"''^. 

However, it is not easy to determine the distribution of aYtaY^ for general 6 E (0, 1): 
indeed let again s„(t) = e'^*T{{aYtaY*)") , n > 1 then 

ri-l 

dtsn{t) = -nY,Sn-k{t)sk{t) + e"*[2n(20 - 1) + {n - l)(n - 2){26 ~ 1)% n>2, 
fe=i 

si{t) = l + (20-l)2(e*-l). 

5.2. Parameter values A G (0, 1],^? = 1/2. In this paragraph, we shall prove the 
following decomposition: 

Proposition 4. Let 

oo 

Mt(z) := ^m„(t)z", ^ > 0, |z| < 1, 

be the moment generating function of the free Jacobi process associated with the pa- 
rameter values X G (0, 1), 9 = 1/2 (we omit the dependence on A for sake of clarity). 
Let also M^o be the moment generating function of the stationary distribution of 
the free Jacobi process associated with the same parameter values (the distribution 
of PUQU*P when U £ is Haar distributed and is -k-free from {P, Q} ). Then 

Mt{z) = M^{z) + ^-^P2Ai/(2-A)[(2 - A)e-*a(z)] + Ut{z) 

where Ut is analytic around the origin whose first and second derivatives vanish, 
and 

lini Ut{z) = 0. 

A->-l 

Proof: Recall that for general parameter values {X,9), the function 

{t,z)^G,{z)='-Y^^ = U('- 
z ^ — ' z" z \ ^ 

n>0 ^ 

defined for t > and |z| > 1, is a solution of the p.d.e. 

dtGt = { [(1 - 2A6i)z - 9{l - X)]Gt + X9z{z - 1)G?} 
with initial value Go{z) = l/(z — 1). Recall also from [9 that 



„ , , (2 - r)z + (1/A - 1) + Vr^z^ -Bz + C' 
^°°^'> 2z(z-l) 

where r = (1/A6'), B = 2{r + {r - 2)/A),C = (1 - 1/A), is the Cauchy-Stieltjes 
transform of the stationary distribution of the free Jacobi process. Then 

{[(1 - 2X9)z - 6{1 - A)]Goo + X6z{z ~ 1)G^} = 
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as can be readily checked from the expression of Goo- In particular, for 6' = 1/2, A G 
(0, 1] the p.d.e 



dtGt = {(1 - A)(2z - l)Gt + \z{z - \)G\) 



holds, while 



•= 2A.(^-1) = -z^' 



satisfies 

{(1 - A)(2z - 1)G^ + Xz{z - 1)G^} = 0. 

Set 

Hf := G( — G(x) 

then 

dtHt = \\z{z - \)El + V422-40 + (l-A)2ift} . 
In particular, if ^ e (0, 1] and 

St[,z) := -i/t f = Mt{z) - M^{z) 



z \z^ 
then 

dtSt = -\d, {a(1 - z)Sl + v'4-4^ + (l-A)2^25i} . 

Now denote 

2 

^*(^) •= ^31^P2At/(2-A) [(2 - A)e a{z)\ 

where we recall that 



Then the p.d.e 



together with 



show that 



dtPt + ^d^pl = 



a'{z) 1 



a{z) z\/\ - z 



dtvt{z) + ^zVT^d,v^{z) = -2e-*a{z){d,pt)[i2 - X)e-'a{z)]. 
Accordingly, the function 

{t,z) ^ ut{z) := St{z) - -j^=vt{z) = ^Ht Q - -j^=Vt{z) 
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satisfies 

dtut = {a(1 - z)ul + 2XVT^utVt + \/4 - 4z + (1 - A)2z2«t+ 

^4 + (l-V^J(.)-^ ' 



(^) 



z 



{a(1 - z)u^t + 2AV1 - zutvt + y/A-Az + il- X)^z^ut+ 



4 + (1 - A)2^.;,| (z) + (a,p,)[(2 - X)e-'a{z)] 

= -^d, {a(1 - z)u^t + 2XVT^utVt + ^4 - 4z + (1 - Xfz^ut+ 
^4+{l-Xr-^vt I (z) + 2ze-'a'{z){dM2 - A)e-*a(^)]. 



Set 



then 



(z) :=^4+(l-A)Y 



(1 - A)2 z2(2 - z) 



^d,{rtvt){z) = - 4^^(^j (l -z)2 ^*^^^ " zr*(z)a'(z)e-*(a,p,)[(2 - A)e-*a(z)]. 
Rearranging terms, one gets 

Stut = -^d, {a(1 - z)u2 + 2VT^UtVt + V4-4^+(l- A)2z2ut} 

+ ^T-tT^T^^*^^) + ^^"'(^tW - 2)a'(z)(5,p,)[(2 - A)e-'a(z)]. 
Wt{Z) {I — z)^ 

Since 

lim rt(z) = 2 

A->-l 

then one can see that 
^* ^= + ^e-*(rt(^) - 2)a'(z)(a.pt)[(2 - A)e-*a(z)] ^ 

as A — )■ 1. As a matter of fact, one can expand 

oo 

Zt = (l-A)^d„(t)z" 
for some coefficients depending on A. Now write 

oo 

= ^^n{t)z^ 
n=l 

and recall the expansion 
Then for A 9^ 1 

00 

V4 - 4z + (1 - A)2z2 = InZ'^ 
n=0 
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where 



7n = 2}_^ 



k=0 ^1^2 



and 

_ 2 + 2v'A(2- A) _ 2-2^A(2- A) 
~ (T^Ap ' ~ (T^Ap ■ 

Besides, 

oo / n \ 
V4-4^+(l-A)222y,(2) = ^K]cfe(i)7n-fe 

n=l \fc=l / 



where we set 



^ 00/ n— 1 \ 

^ n=2 \fe=l / 



^ 00 



n=l 

with 

:= ^ E ^(2 - A)'=-4-i (^fc) e-'^*, n > 1. 

Consequently, 

ciW = -ci{t) + {l-X)di{t). 

But 



1 / (1 - A)(^ - 2) + ^4 - 4z + (1 - A)2z2 \ 1 

Vo(Z) 



l-z \ 2A(l-z) J y/l 

/I 1-A 1 , -, TTT^X 1 1 

( A - - ^V4-4.+ (l-A)2.2j __ _ 

">0 / n>0 n>l 



A 2A 



n=l \ fe=0 / 



In particular 
which implies that 



1 + A 1 1 „ 
^^(°) = ^-2A-2=« 



ci(t) = (1 - A) di(s)e''ds^ e"*. 
Since di (t) = then ci (i) = and 

•00 



n=2 
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Accordingly 

4(0 = -2c2(<) + (1 - X)d2{t), 

and since 



then 



But after some effort, one sees that d2{t) — yielding C2{t) — and 



oo 



71 — 3 

Similarly 

= -3c3(t) + (l-A)d3(i) 

and 

(1-A)(3-A) 
= 32 • 

It follows that 

C3(t) = (1 - A) (^-^ + ^* d3(s)e3^ds) 
and computer assisted computations show actually that 

C3W = -^[2Ae-3* + 3(l-A)e-*] . 

More generally 

Cn n-2 
X!cfcW7n-fc + 2V'i(t)c„_i(t) + ^ [A(c„_fc(t) - C„_i_fcW) 
fc=3 fc=3 

(11) +2(V„-fe(i) - ^«-i-fe(i))] cfc(i)) + (1 - A)d„(0. 

for n > 4, where the last sum is taken to be zero when n — \. Finally, it remains 
to prove that 

lim c„(t) = 0. 

A-i-l 

To proceed, it is easy to note from (|10p that 

c„(0) =c„_i(0) + fc„(A) 

for any n > 2, where 
fc„(A) = V«_i(0)-Vn(0)-^7« 

ri-\ 



?2ri-l ' 22 



^^^Y4/^ 2n-2 \ / 2n \\ 1 ^ /3fc/3„_fc 



k=l V " 'V V" 'V/ '"fc=0^1^2 



But 



together 



0^ 1 _ (1-A)^^(1 + VA(2^)"-^^ 

— ^k^n-k on — 



hni V4-4z+(l-A)2z2 = 2^/l~7 

A->-l 
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entail 

hPn-k _ -2(2n)! 



iim 7„ — lim z > — r — , ^ , ^ . 

A^i A^i ^ zfz?"'^' 2n-l (2"n!)2 

1 1 ^V./^ 2n-2 \ / 2n \\ (2n) 



fe=o -'I -'a 

As a result 

n-l 



A^i ' 22»-i 22»-i V V'^-l-V \n-kJJ (2?i- l)(2«n!)2 

and applying twice Pascal rule shows that 

2n \ / 2n-2 \ /2n-2\ 2n-2 
— k J \n — k — 2 ) \ n — k J \n — k — 1 

Hence 

2n-2 \ / 2n \ /2n-2 \ /2n-2 \ /2n-2 \ /2n-2 



yn— 1 — fcy \n — k J \n — k — 1) \n ~ k — 2 J \n — k — 1 J \n — k 
yielding 

, , 1 ( (2n-2\ (2n-2\\ (2n)! 

As a matter fact, it suffices to use induction in order to prove that c„(0) — > 0, n > 4 
as A — > 1. Once this is proved, it is clear from (fTTj) that 



lim c„(f) = 

A->-l 



for any n > 4 and any t > 0. 



5.3. Parameter values A G (1,2),0 = 1/2. As mentioned in the introduction, 
the moments of the free Jacobi process associated with the parameter values A G 
(1,2),6' = 1/2 may be deduced from those corresponding to A G (1,2),6' = 1/2. 
Indeed 

Proposition 5. 



r[((l - P)YtQY:T\ = r{Q) + ^(-1)'^ UiPY.QYrf)- 

k=l ^ ^ 

Proof. Let n> \ and define for oi, 02, . . . , a„ G 

/i„(ai, . . . ,a„,i) := T{aiYtQY^ . ..a^YtQY^). 
hn is obviously a linear map with respect to its n first coordinates, thereby 
hn{l - P, ■ . ■ ,1 - P,t) ^ ^ /i„(ai^, • • • ,aj„,i) 

(Qi,,...,ai„)6{l,-P}" 

Since Q is idempotent and r is a trace, then one realizes that each n-tuples with 
k,l < k < n elements equal to (—P) and n — k elements equal to 1 contributes to 
{—l)'^T{{PYtQYf')^). They are even the only tuples giving this contribution since 
there is no relation between P and YtQY*. Since their total contribution is 

i-l)' (^^T[{PYtQY;)% l<k<n. 

while that of the 7i-tuple (1, . . . , 1) is t(Q), then the proposition is proved. I 
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